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We present bulk property measurements of NpIr, a newly synthesized member of the Np-Ir binary 
phase diagram, which is isostructural to the non-centrosymmetric pressure-induced ferromagnetic 
superconductor Ulr. Magnetic susceptibility, electronic transport properties at ambient and high 
pressure, and heat capacity measurements have been performed for temperatures T = 0.55 — 300 K, 
in a range of magnetic helds up to 14 T and under pressure up to 17.3 GPa. These reveal that Npir 
is a moderately heavy fermion Kondo system with strong antiferromagnetic interactions, but there is 
no evidence of any phase transition down to 0.55 K or at the highest pressure achieved. Experimental 
results are compared with ab initio calculations of the electronic band structure and lattice heat 
capacity. An extremely low lattice thermal conductivity is predicted for Npir at temperatures above 


300 K. 

PACS numbers: 61.05.cp, 72.15.-v, 75.20.Hr, 75.40.Cx 

I. INTRODUCTION 

The magnetic and electrical properties of the major¬ 
ity of rare earth intermetallics are well explained using 
the standard localised moment model of rare earth mag¬ 
netism, as expounded by Jensen and Mackintosh^, ow¬ 
ing to the limited radial extent of the 4/ wave-functions. 
Whereas, upon descending down the periodic table into 
the 5/ actinide series, the /-electron wave-function be¬ 
comes more extended. Therefore, the 5/ electrons have 
a character intermediate between that of the localised 
4/ electrons and the itinerant 3d electrons of the tran¬ 
sition metals, which taken in conjunction with the in¬ 
creased difficulties in handling such materials, results in 
a more limited understanding of magnetism in the ac¬ 
tinides. However, this also contributes to the fact that 
actinide systems display a unique complexity, exhibiting 
various interesting properties such as heavy fermion and 
non-Eermi liquid behaviours^, multipolar order— and 
unconventional superconductivity^^— . 

Ulr has recently been identified as a ferro¬ 
magnetic quantum critical point pressure-induced 
superconductor— and belongs to two different, non- 
conventional subclasses of superconductor: non- 

centrosymmetric superconductors such as CePtaSii^, and 
ferromagnetic superconductors such as UGe 2 ^. It is 
thought that the absence of inversion symmetry prohibits 
spin-triplet (p-wave) pairing, whilst ferromagnetism pro¬ 
hibits spin-singiet (s-wave) pairing. Therefore Ulr is a 
particularly interesting system, significant for the study 
of the interplay between magnetism and superconductiv¬ 
ity in strongly correlated electron systems. It is an itiner¬ 


ant ferromagnet (Tq = 46 K at ambient pressure), with 
an ordered moment of only O.G/rs/U atontii. However, 
at high temperatures its behaviour is better described in 
a localised picture, with an effective moment of 3.57/iB 
(consistent with either 5/^ or 5/®) obtained from a Gurie- 
Weiss fit to single crystal high temperature magnetic sus¬ 
ceptibility measurements (T = 300 — 800 K}^. 

Studies of isostructural transuranium compounds offer 
the opportunity to investigate how physical properties 
evolve as a function of 5/ shell filling, with the possibil¬ 
ity in the case of Npir of investigating the significance of 
the itinerant ferromagnetism to the nature of the super¬ 
conductivity. Here we report on the synthesis and char¬ 
acterisation of Npir, an isostructural analogue of Ulr, 
via magnetic susceptibility, resistivity and heat capac¬ 
ity measurements. The experiments reveal no magnetic 
ordering or superconductivity down to 0.55 K, in con¬ 
trast to the development of ferromagnetic order in Ulr 
at Tc = 46 K, and no evidence of pressure-induced su¬ 
perconductivity. Ab initio calculations reveal the signifi¬ 
cance of the Hubbard -PU term for the description of the 
Np 5/ electrons. The calculated phonon dispersions, lat¬ 
tice heat capacity, and lattice thermal conductivity are 
presented. The latter is found to be exceptionally low in 
the high temperature regime. 


II. METHODS 

Poly crystalline samples of Npir were synthesised at 
the Institute for Transuranium Elements (ITU) by arc 
melting stoichiometric amounts of Np (99.9%) and Ir 
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FIG. 1: (Color online) Rietveld refinement (solid red line) 
of room temperature x-ray powder diffraction data for Npir 
(o) annealed at 873 K. The difference between the calculated 
and experimental points is shown by the black line which 
has been offset by -2000 cts/4s for clarity. The vertical tick 
marks correspond to the Bragg peak positions for the NpIr 
P2i structure shown in the inset. 


(99.98%) metals under a high purity argon atmosphere 
(Ar: 6N) in a water cooled copper hearth using a zir¬ 
conium getter. The sample was re-melted five times to 
obtain a good homogeneity and the weight loss was 0.5%. 
The as-cast sample was embedded in a tantalum foil and 
encapsulated in a quartz tube under high vacuum, and 
then annealed at 873 K for 3 weeks. X-ray powder diffrac¬ 
tion analysis of the samples was performed on a D8 Ad¬ 
vance diffractometer with a Bragg-Brentano configura¬ 
tion, equipped with a Cu X-ray tube (40 kV, 40 mA), a 
Ge monochromator (111), and a Lynx Eye linear position 
sensitive detector. The powder pattern was recorded at 
room temperature in step scan mode over a 29 range of 
10 — 120°, with a step size of 0.013° and a count time of 
4 s per step. Npir was indeed found to be isostructural to 
Ulr, the x-ray diffraction pattern being fitted with a mon¬ 
oclinic structure described by the P2i space group, with 
lattice parameters a = 5.5832(6) A, b = 10.7368(9) A, 
c = 5.5848(6) A and /3 = 95.708(5)° (see Table U for re¬ 
fined atomic parameters). Figure [T] shows the quality of 
the refinement. At the lowest 26 angles there is a broad 
peak arising from the necessary encapsulation of the sam¬ 
ple, but the other features of the data are reproduced by 
the refined structure. 

Magnetic susceptibility and isothermal magnetisation 
measurements were performed using a Quantum Design 
MPMS-7 Squid magnetometer on a 89.9 mg Npir sample 
over a temperature range of 2 — 300 K, in magnetic fields 
up to 7 T. 

The ambient pressure electrical resistivity, magnetore¬ 
sistivity and Hall effect have been measured in the tem¬ 
perature range 1.8 — 300 K, and in magnetic fields up 
to 14 T, using a Quantum Design PPMS-14T setup, by 
means of a 4 DC probe technique voltage measurement. 


TABLE I: Refined atomic parameters for Npir (space group 
P2i) at room temperature, with i?„p = 10.52 and GoF = 
2.14. 


atom 

Wyckoff 

X 

y 

z 

Occ 

Npl 

2a 

0.13759 

0.00000 

0.12635 

1 

Np2 

2a 

0.62272 

-0.00257 

0.63025 

1 

Np3 

2a 

0.87689 

0.71522 

0.38225 

1 

Np4 

2a 

0.39253 

0.71821 

0.87415 

1 

Irl 

2a 

0.10633 

0.26729 

0.10832 

1 

Ir2 

2a 

0.62036 

0.25888 

0.62317 

1 

Ir3 

2a 

-0.12913 

0.44782 

0.37114 

1 

Ir4 

2a 

0.35812 

0.45145 

0.85769 

1 


An Npir sample of size ~ 1.5 x 0.4 x 0.2 mm^ was pol¬ 
ished on two parallel faces to determine better the form 
factor. Electrical contacts between the sample surfaces 
and the 50 /rm silver wires were ensured by using silver 
epoxy (Dupont 4929). Finally, each mounted sample was 
then encapsulated with Stycast epoxy (1266). For electri¬ 
cal resistivity measurements, the current I was applied in 
the polished plane. For the magnetoresistivity, I was par¬ 
allel to the voltage direction and parallel to the applied 
magnetic field B. In the Hall configuration, the voltage 
Vh was measured perpendicular to the current / and the 
applied magnetic field B. The Hall resistance (Rh) has 
been determined by measuring Vh under fields alternat¬ 
ing between -1-14 and —14 T. The magnetic field response 
Vh{B) at fixed temperatures has been measured to con¬ 
firm results obtained when ramping in temperature. For 
all measurements 7 = 5 mA was used. 

The high-pressure resistance measurements were per¬ 
formed by a four-probe DC method in a Bridgman-type 
clamped pressure cell, with a solid pressure-transmitting 
medium (steatite). Electrical contacts were made with 
25 p,m diameter platinum wires lightly pressed onto the 
sample. Before each measurement the cell was loaded and 
clamped at room temperature. The exact pressure inside 
the pressure cell was determined later by using the pres¬ 
sure dependence of the superconducting transition tem¬ 
perature of lead as a manometer—. Measurements were 
performed on a ^ 50 fig sample of size ~ 20 x 50 x 500 fim^ 
in the temperature range 1.8 — 300 K up to 17.3 GPa. 

Heat capacity measurements on a 5.6 mg sample of 
Npir were performed over the temperature range T = 
2 — 250 K, and on a 0.9 mg sample down to 0.55 K, 
via the standard relaxation calorimetry method using a 
Quantum Design PPMS-9 with the ^He refrigeration in¬ 
sert, after the samples had been coated in Stycast. The 
data have been corrected for the addenda and the sty¬ 
cast. Self-heating effects in neptunium make it difficult 
to reach lower temperatures. No isostructural phonon 
blank exists, so instead, in order to best estimate the 
phonon contribution, we have followed two strategies. 
First, we have synthesized orthorhombic Thir (space 
group Cmcm) and cubic Luir (Pm3m), verified their 
structures and phase purities using x-ray powder diffrac¬ 
tion, and measured their heat capacities using the same 
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relaxation method. Second, we have calculated ab initio 
the phonon spectrum of Npir and the lattice contribution 
to the heat capacity. 

III. RESULTS AND DISCUSSION 
A. Magnetic Susceptibility 

The magnetic susceptibility data for NpIr, shown in 
Fig. Ufa), reveal no indication of any magnetic transi¬ 
tion above T = 2 K. There is no difference in the re¬ 
sults for zero-field cooling and in-field cooling. Between 
50 and 300 K the inverse susceptibility may be mod¬ 
elled by a Curie-Weiss law, see the inset to Fig. Elja), to 
obtain a Curie-Weiss temperature of —104 ± 1 K, and 
an effective paramagnetic moment of 3.01 ± 0.02 /i^/Np. 
Starting with the Curie-Weiss temperature, such a large 
negative value is indicative of strong antiferromagnetic 
interactions. Intriguingly, this is of a similar magnitude 
to the —300 K obtained from high temperature measure¬ 
ments on Ulr— , but in that case, in spite of the anti- 




hoH (T) 

FIG. 2: (Color online) (a) Magnetic susceptibility of Npir 
measured with /roR = 1 T. The inset shows a Curie-Weiss fit 
to the inverse susceptibility, (b) Isothermal magnetisation of 
Npir measured at T = 2,10, and 40 K. 


ferromagnetic interactions, it undergoes a ferromagnetic 
phase transition, highlighting that these compounds sit 
at the interface between antiferromagnetic and ferromag¬ 
netic order. Next we consider the value for the effective 
paramagnetic moment, which is inconsistent with either 
a 5/^ or a 5/^ configuration in both the Russell-Saunders 
(3.62 /xb/Np and 2.68/iB/Np) and the intermediate cou¬ 
pling schemes (3.68/iB/Np and 2.76 ^s/Np^S), and sug¬ 
gests, therefore, that possibly the 5/ electrons are not 
fully localised. If alternatively the fit to the inverse sus¬ 
ceptibility is only made above T = 250 K, then a Curie- 
Weiss temperature of —155 ±5 K and an effective param¬ 
agnetic moment of 3.20 ± 0.02 /r^/Np are obtained, still 
inconsistent with a localised moment picture. Such re¬ 
sults are not wholly dissimilar to those for Ulr, for which 
a localised behaviour was only observed at high temper¬ 
atures (r > 300 K)i^. Regrettably, it was not possi¬ 
ble to measure the magnetic susceptibility above room 
temperature using our experimental setup. However, as 
might be expected given the standard trend towards in¬ 
creasingly localised electrons on spanning the actinide 
series, Npir is perhaps more localised at low tempera¬ 
tures than Ulr, as demonstrated by the difference in the 
effective moments, obtained by making a Curie-Weiss fit 
to the inverse susceptibility over the temperature range 
T = 50 - 100 K, of 1.67^B/Ui2 and 2.96 ± O.Ol/rs/Np 
for polycrystalline Ulr and Npir, respectively. 

Fig. [2Kb) presents the isothermal magnetisation data 
for T = 2,10 and 40 K. These reveal no saturation for 
magnetic fields up to 7 T, and no evidence of magnetic 
hysteresis. 


B. Resistivity 

The electrical resistivity of an annealed polycrystalline 
sample of Npir is shown in Fig |3Ka). Interestingly, 
the room temperature absolute resistivity of Npir is 
p = 122 /iHcm, which is comparable to the 80 ^Hcm 
reported for the best quality Ulr crystals^!. With de¬ 
creasing temperature, we observe the presence of a very 
broad maximum centred at T„iax 150 K, reminis¬ 
cent of that observed in several neptunium-based sys¬ 
tems such as NpCoGag^ and NpPdj^, which is indica¬ 
tive of a Kondo-type behaviour. Below 50 K coherence 
sets in and the resistivity collapses to only 1.95/iHcm 
at T = 1.8 K. Below 6 K (see inset to Fig. |3Ka)), the 
strong curvature can be modelled by a Fermi liquid be¬ 
haviour: p = po + with po = 1-42 ± 0 . 02 /iHcm 

and Afl = 0.209 ± 0.001 pHcmK“^. This gives a resid¬ 
ual resistivity ratio prt/pq ~ 90, which approaches the 
value of 230 reported for highest quality single crystal 
Ulr—. Over the temperature range T = 5 — 20 K, 
the data can be modelled by a non-Fermi liquid law: 
p = po + with po = 0.92 ± 0.07/iHcm and 

A 5/3 = 0.412 ± 0.001 pOcmK“^. Such a variation in the 
temperature dependence with temperature interval may 
indicate different spin fluctuation regimes within Npir. 
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FIG. 3: (Color online) (a) Electrical resistivity of Npir, with inset showing fits to a Fermi liquid (solid red line) and non-Fermi 
liquid (dashed green line) model at low temperatures, (b) Isofield magnetoresistance of NpIr at 14 T, with inset showing the 
field dependence of the isothermal magnetoresistance for T = 1.8 K (o), 3 K (□), 5 K (A), and 50 K (*). (c) Hall Effect 
measurements of Npir, with inset showing fit to Rh{T) = Ro + Rix*{T). 


Fig. |3Kb) shows the longitudinal magnetoresistivity of 
Npir at 14 T over the temperature range T = 1.8 — 50 K, 
which is similar in shape to that of UAI 2 below 20 K— . 
For all temperatures up to 50 K the magnetoresistive 
contribution is positive and, at T = 1.8 K, large relative 
to the resistivity (~ 140%). There is no clear evidence 
of any anomaly down to this temperature that might be 
associated with a magnetic phase transition. The field 
dependence of the isothermal magnetoresistivity for a 
range of different temperatures is shown in the inset to 
Fig. 13(b). For all temperatures measured (T = 1.8, 3,5 
and 50 K), the magnetoresistivity varied quadratically as 
a function of the magnetic field. 

Figure EKc) presents the Hall coefficient {Rh) for T = 
1.8 — 75 K for Npir. Rh is slightly enhanced approach¬ 
ing 10 X 10“^°m^C“^ at 50 K, which is nevertheless 
rather similar to values for metallic systems, e.g. Cu: 
0.5 X 10“^°m^C“^. Upon decreasing the temperature, 
Rh decreases almost linearly, changing sign at 25 K, be¬ 
fore reaching a minimum at 7 K, below which Rh starts 
to increase again. As Rh is the result of the combina¬ 
tion of the electron and hole contributions with different 
carrier velocities and relaxation times, our data suggest 
that the nature and the mobility of the carriers are chang¬ 
ing drastically with temperature. The Hall coefficient is 
composed of two terms: Rh = Ro + Rs, where Rq is the 
ordinary Hall coefficient, and Rs is the extraordinary or 
anomalous Hall coefficient. Following Ref. we replace 
Rs by a term dependent on the magnetic susceptibility 
to give: 

RHiT) = Ro+Rix*iT), ( 1 ) 

where the first term, Rq, describes the Hall effect due to 
the Lorentz motion of the carriers and/or residual skew 
scattering by defects and impurities, while the second 
term comes from skew scattering by Hondo impurities. 
In this formula x* is the reduced susceptibility, approxi¬ 
mated by x(T)/C, where C is the Curie-Weiss constant 
obtained from the fits to the inverse susceptibility above. 
Assuming Rq and Ri are independent of the temperature. 
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FIG. 4: (Golor online) (a) Electrical resistance of Npir mea¬ 
sured as function of temperature for a selection of applied 
pressures, (b) Comparison of normalised residual resistance 
data for Npir for pressures between 0.16 GPa and 17.65 GPa. 
(c) Temperature exponent extracted from low temperature 
fits to the resistance of 7? = AT" -|- Rq. 


plotting Rh{T) as a function of X*(T), as shown in the 
inset to Fig.[3Kc), gives Rq = -|-7.19±0.02 x 10“® m^C“^, 
indicating that the ordinary Hall effect is dominated by 
the hole contribution. A simple one-band model then 
provides an estimation of 8.69 ± 0.03 x 10^® m“^ for the 
concentration of free holes, giving an upper limit for the 
actual carrier concentration in Npir in the normal state. 
This may then be converted into a rough estimate of 0.05 
for the number of free holes per formula unit at high tem¬ 
perature. 

The global shape of the electrical resistivity of Npir 
does not change drastically under pressure up to 12 GPa 
(Fig. |4|). However, with increasing pressures up to 
^12 GPa, the maximum around 150 K, seen in Fig. [3(a), 
becomes less apparent, but the onset of coherence below 
this temperature remains. Above ~12 GPa, the resis- 
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tivity evolves more smoothly with temperature with a 
shift of the scattering to higher temperature. At low 
temperature the resistance can be fitted according to 
P = Po + alT", where the exponent stays essentially con¬ 
stant ^1.33±0.15 until above 13 GPa, where it fluctuates, 
as shown in Figure lUJc). No hint of a superconducting 
transition is detected down to 1.8 K for any pressure be¬ 
low 17 GPa. 

One possible explanation for the observed absence of 
ferromagnetism in Npir above 1.8 K may be inferred from 
a comparison of the forms of the resistance curves of 
Ulr and NpIr. At ambient pressure, the resistivity of 
NpIr (Fig. [SJa)) resembles the resistivity of Ulr under 
pressures greater than 2 GPai^i, which appear to be un¬ 
favourable conditions for ferromagnetism in Ulr. The 
requirements for the appearance of superconductivity in 
Ulr are extremely drastic, and strongly dependent on the 
pressure transmitting medium^ij but our measurements 
on Npir were performed under similar hydrostatic con¬ 
ditions. The superconducting transition temperature of 
Ulr is lower than our accessible temperature range, so it 
is possible that Npir may still superconduct at tempera¬ 
tures below 1.8 K, but the lack of a ferromagnetic state, 
which seems to be a prerequisite for non-conventional su¬ 
perconductivity in Ulr potentially makes this less likely. 


C. Heat capacity 

As shown in Fig. [5l the heat capacity of Npir varies 
smoothly between 2 K and 270 K, with no anomalies 
which might be associated with any phase transition. 
When a straight line fit is made to Cp/T versus 
for T <7 K, as in the inset to Fig. [SJ we obtain 7 = 


60 
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FIG. 5: (Color online) Heat capacity of Npir, Thir and Luir, 
with inset showing a fit to CpjT = 7 -I- /3T^ at low tempera¬ 
tures. Also shown as a dashed line is the Dulong-Petit limit of 
QR for the phonon heat capacity, and the ab initio calculated 
lattice specific heat as a solid black line. 



T(K) 

FIG. 6 : (a) Total electronic (5/-l-conduction electrons) con¬ 
tribution to the heat capacity of Npir divided by temperature 
obtained from + 7 '^*'^" • T. (b) The entropy ob¬ 

tained by integrating the data in panel (a) as a function of 
temperature. 


175 ±1 mJK“^mol“^ for the electronic heat capacity and 
a Debye temperature Qp = 145.7±0.5 K. The electronic 
heat capacity of Npir is considerably greater than that 
for Ulr (7 = 49 mJK“^mol“^, Ref. [2^, and is indicative 
of strong electronic correlations. Gombining the value 
of 7 with the coefficient Afl for the quadratic term in 
the low temperature resistivity obtained in section IIIIBl 
gives a Kadowaki-Woods ratio^l Afl/t^ = 0.68 ± 0.01 x 
10 “^ /rDcmK^mol^mJ”^, implying that Npir is a mod¬ 
erately heavy fermion material, which would be consis¬ 
tent with hybridisation causing an effective paramagnetic 
moment below the free ion value. Such a value for the 
Kadowaki-Woods ratio is comparable with that for UAI 2 
(0.89 X10“^ ^OcmK^moPmJ“^)21, in which the low tem¬ 
perature specific heat is well expressed in terms of spin 
fluctuations, that prevent any magnetic order above I K. 
However, the paramagnon upturn present in UAI 2 data is 
absent in that for Npir. USns also displays a similar value 
of Afl/7^ = 0.78 X lO-^GcmK^moPmJ-^^ with a 
very similar specific heat value, 7 = 172 mJK“^mol~^ 
(Ref. [ 28 I) and despite the lack of a paramagnon upturn 
in the low temperature heat capacity, is also classified as 
a spin fluctuator system. 
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Figure [5] also compares the heat capacity of Npir with 
that of Luir and Thir, in an attempt to estimate the 
lattice contribution. Regrettably, neither ThIr nor LuIr, 
crystallising in the Cmcm and PmZm spacegroups re¬ 
spectively, are isostructural with NpIr. However, we have 
chosen to use the data for Thir as a phonon blank, since 
its molar mass is very similar to that of Npir, and its De¬ 
bye temperature (Qd = 176.6±0.3K) is closer than that 
of Luir [Qd = 254.7 ±0.4 K). Furthermore it agrees bet¬ 
ter with the ab initio calculated (see Section ITVT) phonon 
heat capacity of Npir (shown as the black solid line). 
Hence we estimate the total electronic (5/ ± conduction 
electron) contribution to the heat capacity of Npir as: 

Cp = — Cp^^"^ (phonons), 

Cf = C'pP^"'- ± 7™’'• T, (2) 

where 7 ^^’' = 4.8±0.1 mJK“^mol”^ is the electronic heat 
capacity of Thir. Fig. [SKa) displays the total electronic 
contribution to the heat capacity of Npir divided by tem¬ 
perature, revealing a broad peak centred at T = 30 K, 
and a near constant behaviour at high temperatures. 

Integrating the total electronic contribution to the heat 
capacity of Npir allows an estimate for the entropy to 
be obtained, which is shown in Fig. [61(b). The entropy 
varies smoothly, and by extrapolation to higher tempera¬ 
ture appears to be compatible with the free ion value for 
Np^+ (RlnlO) or Np"‘+ (i?ln9). Assuming that the 5/ 
contribution to the heat capacity is close to saturation by 
~ 200 K, we can deduce from Fig.lHKa) that the electronic 
coefficient has an upper limit of 7 = 20±2 mJK“^mol“^ 
at high temperature. 

Such a considerable difference between the low and 
high temperature electronic heat capacities suggests that 
the enhanced low temperature 7 value is due to a strong 
Kondo interaction while localisation features are en¬ 
hanced at higher temperatures. 

IV. AB INITIO CALCULATIONS 
A. Methodology 

Electronic structure calculations were carried out using 
the Vienna Ab-initio Simulation Package (VASP) ^^’^° , 
with the generalized gradient approximation (GGA) 
as the density-functional theory (DFT) exchange- 
correlation functional, as well as with its extension to 
treat strongly correlated electrons, DFT with an ad¬ 
ditional Hubbard U term (DFT±17) ^^i^^ . Within the 
GGA±C/ approach, we have used the Dudarev et al. 
formulatioii^, where the Hubbard and exchange param¬ 
eters, U and J, respectively, are introduced to account 
for the strong on-site correlations between the neptunium 
5/ electrons. This helps to remove the self-interaction er¬ 
ror and improves the description of correlation effects in 
the open 5/ shell. We have chosen a Hubbard U value 
of 4.0 eV and an exchange parameter J value of 0.6 eV, 


FM PM AFM 



FIG. 7: (Color online) Atom-resolved partial density of states 
for the three different magnetic orders calculated for Npir 
using the GGA-|-[/ approach. 


which are in the range of accepted values for Np and Pu 
compounds^i^. To test the dependence of our results 
on the U value, we have also performed calculations for 
U = 2 and 3 eV, and for 1/ = 0 eV, i.e., for the com¬ 
mon GGA functional. Further, to deal with the problem 
of degenerate metastable states when using the DFT±C/ 
methodology, we have used the occupation matrix control 
(OMG) method proposed by Dorado et alM- This method 
consists of the direct control of the strongly correlated 
electron occupation matrices. Details of the electronic 
structure calculations can be found in Ref. [s^ 

We have considered three different magnetic or¬ 
ders: ferromagnetic (FM), antiferromagnetic (AFM), 
and paramagnetic (PM) order. In the FM ordered state, 
we assume that all the Np ions have collinear mag¬ 
netic moments oriented along the c-direction. In the 
AFM ordered state, the Np ions are considered to be 
collinear with magnetic moments changing sign from 
one Np plane to another. Finally, for the PM ordered 
state, we adopt the disordered local moments (DLM) 
approac h^^i^^ , which states that paramagnetism can be 
modelled as a state where atomic magnetic moments 
are randomly oriented (noncollinear magnetism), valid 
for materials that display a Gurie-Weiss paramagnetism, 
such as Npir. The DLM approach can be simplified by 
considering only collinear magnetic moments when the 
spin-orbit coupling is not taken into account. Hence, the 
problem of modelling paramagnetism becomes a prob¬ 
lem of modelling random distributions of collinear spin 
components. It can be solved by using special quasi¬ 
random structures (SQS)^. An SQS is a specially de¬ 
signed supercell built of ideal lattice sites to mimic the 
most relevant pair and multisite correlation functions of 
a completely disordered phase (PM order in our case). 
As a PM simulation cell we used an extended lattice cell 
of 64 atoms. We note however that on account of the 
large simulation cell needed for the DLM calculations, it 
was not possible to perform these including the spin-orbit 
interaction. 
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B. GGA+{7 results 

A full structural and atomic-site relaxation has been 
carried out for Npir (with [/ = 4 eV and J = 0.6 eV). We 
have found, in agreement with experiments, that NpIr 
crystallizes in a monoclinic structure (described by the 
P2i space group). We further found that the calcu¬ 
lated lattice parameters a = 5.6072 A, 6 = 10.7829 A, 
c = 5.6088 A and jS = 95.708° are very close to the ex¬ 
perimental values (see Section P. The relaxed atomic 
positions for the FM order are given in the Appendix „ 
where again the proximity to the experimental atomic 
positions can be observed. 

The total energy for the three different magnetic orders 
have been calculated. Although the energy differences 
are small, we found that the FM order has the lowest 
total energy, followed by the PM and AFM orders, having 
0.039 and 0.056 eV per formula unit higher total energies, 
respectively. These findings differ from the experimental 
results which reveal no sign of magnetic order for T > 
1.8 K, while the Curie-Weiss temperature implies strong 
antiferromagnetic interactions. 

To investigate the origin of the obtained energy or¬ 
der we have investigated the influence of the Hubbard U 
parameter and considered the influence of the spin-orbit 
interaction. Performing calculations with U values of 0, 
2, and 3 eV did not lead to a change in the relative en¬ 
ergy sequence. The FM phase was always found to have 
the lowest total energy, followed by the PM phase, and 
then by the AFM phase. The relative energy differences 
were similar to those found for [/ = 4 eV. We emphasize, 
however, that these calculations were performed without 
the spin-orbit interaction, as the DLM calculations are 
computationally too heavy with spin-orbit interaction. 
We have for sake of comparison computed a hypotheti¬ 
cal nonmagnetic phase of Npir (i.e., no moments at all) 
with the spin-orbit interaction. We find its total energy 
to be higher than that of both the FM and AFM phases. 
This indicates that a nonmagnetic state is unlikely and 
would be an insufficient representation of the PM phase. 
Thus, in the absence of DLM calculations with spin-orbit 
interaction we cannot definitely state what the lowest en¬ 
ergy magnetic order of Npir is It is however worth noting 
at this point that isostructural Ulr, which also displays 
antiferromagnetic interactions based on the Curie-Weiss 
temperature, in fact orders ferromagnetically below 46 K. 
The local spin moments of the Np ions have also been 
calculated; these are 3.78 /tb, 3.74 and 3.71 ^b for the 
FM, PM, and AFM orders, respectively. 

The calculated partial densities of states (DOS) for 
the three different magnetic orders are given in Fig. [T] 
There are significant hybridizations of the Np-/ and Ir-d 
electrons in the energy range of —3 to -1-3 eV, as can be 
inferred from the similar DOS structures. The Coulomb 
U potential leads to a splitting in the 5/ spectrum (of 
about 3 eV) which can be clearly seen for the PM and 
AFM ordered states. 

To assess the importance of the spin-orbit interaction 



-6 -4 -2 0 2 4 


Energy (eV) 


1.2 



-6 -4 -2 0 2 4 

Energy (eV) 

FIG. 8: (Color online) Top: Total and partial density of states 
of Npir calculated using the GGA-tC approach and including 
the spin-orbit coupling. Bottom: The same, but computed 
with the GGA approach and including the spin-orbit coupling. 


on the atomic magnetic moments we have carried out an 
analysis of its influence on the FM order. We find that 
the total magnetic moment M, written as sum of the spin 
magnetic moment Ms and the orbital magnetic moment 
Mq is M = Ms Mq = 3.748/iB — 0.933 ^b = 2.815/iB- 
This value is very close to the intermediate coupling value 
for a 5/^ configuration. The total and partial density of 
states including spin-orbit coupling are shown in Fig. [5] 
(top). 


G. GGA results 

To investigate the importance of the Hubbard U term, 
we have performed a similar study using the plain GGA 
approach. Similarly to the spin-orbit case, we only an¬ 
alyzed the FM order. For this approximation the lat¬ 
tice parameters are a = 5.4620 A, 6 = 10.5037 A, 
c = 5.4636 A and j3 = 95.708°. Although they are in 
good agreement with the experimental results, the devi¬ 
ation with respect to these is larger than when comparing 
with the GGA-|-[/ results. The local magnetic moments 
on the Np ions without accounting for the SO coupling 
have a magnitude of 3.113 /tb. If we include the SO cou¬ 
pling the spin magnetic moment drops to 2.842 /ib while 










the orbital magnetic moment becomes —2.527/is, giving 
a net moment of only 0.315 /tb- The calculated density of 
states are plotted in Figure [8] (bottom), where it can be 
observed that there is a splitting of the / orbitals (into 
5/5 and 5 / 7 ) in the case that SO interaction is included. 
Nonetheless the manifold of 5/ states appears close to the 
Fermi energy and has its maximum at the Fermi energy. 
Application of the GGA+[/ method conversely splits the 
5/ manifold of states and leads to a low 5/ DOS near 
the Fermi level (Fig. HI top). 



Phonon DOS/atom 


D. Phonons properties and thermal condnctivity 

We have calculated the phonons of Npir using the 
finite-displacements method in conjunction with super¬ 
cells consisting of 2 x 2 x 2 primitive cells (128 atoms). 
The interatomic forces were calculated with VASP, 
adopting as above, the GGA-bG approach for the elec¬ 
tronic structure. The phonon modes were obtained with 
the phonopy package^ in the quasiharmonic approxima¬ 
tion. To enable this approximation the system volume 
has been isotropically expanded by 2 % from the GGA-|-t/ 
relaxed volume. The anharmonic effects induced by the 
volume dependence of phonons frequencies are explored 
and the lattice thermal properties such as the lattice spe¬ 
cific heat and the phonon thermal conductivity are cal¬ 
culated. 

The ab initio calculated phonon density of states and 
phonon dispersion curves a;„q are given in Fig. HI The 
16-atom unit cell of NpIr results in 48 phonon modes 
with a rather homogeneous spreading of the bands from 
0 to 4 TFIz. The atom-projected phonon DOS (right- 
hand panel) shows that at low energies the contribution 
from the Ir atoms is larger than that from the Np atoms, 
while at higher vibrational frequencies this behavior is 
reversed. 

The lattice heat capacity Cp can be obtained from the 
Gibbs free energy G{T,p) at constant pressure, Cp = 
—T{d^G/dT^). The Gibbs free energy is obtained from 

GiT,p) = mm [U{V) + V) + pV] , (3) 

where U{V) is the volume-dependent electronic total en¬ 
ergy and the phonon free energy. 


FIG. 9: (Color online) GGA-I-G calculated phonon disper¬ 
sions (left panel) and corresponding projected phonon density 
of states per atom (right panel) of Npir. The q-point labels 
in the left panel are those for the standard high-symmetry 
positions of the monoclinic primitive Brillouin zone. 


Lastly, we investigate the lattice thermal conductivity 
kl of Npir. A low thermal conductivity of materials is of 
interest as this would lead to a high thermoelectric figure 
of merit^. The thermal conductivity of Npir is unknown, 
but, its anisotropic, low-symmetry crystal structure sug¬ 
gest that its lattice contribution could be very small. We 
have computed the (direction averaged) hl{T) of Npir 
using an approximate solution to the phonon Boltzmann 
transport equation in the relaxation-time approximation, 

«:L(T) = i^ f vl^TnqCnci, (5) 


where the sum is over all phonon modes, Vnq is the group 
velocity of a given phonon mode, C„q is the mode heat 
capacity depending only on the mode frequency ujnq and 
the temperature, and T„q is the mode dependent relax¬ 
ation time, which is computed here on the basis of the 
model of Bjerg et al^ Furthermore, for the determina¬ 
tion of the lattice thermal conductivity the Griineisen 
parameter 7 is a fundamental quantity. It characterizes 
the relation between phonon frequency and crystal vol¬ 
ume change, and is defined as 


Vuc duJnq 

^nq 


( 6 ) 


poo 

FPhor^{T,V) = / dujg{oj,V)[haj/2 

Jo 

(4) 

with g{u}, V) the phonon DOS, computed as mentioned 
above for different volumes. The ab initio calculated 
Cp{T) of Npir is shown in Fig. [S] As can be noted the 
computed lattice heat capacity is smaller than the mea¬ 
sured heat capacity of Npir. Its temperature dependence 
corresponds very well with the measured heat capacity 
of Thir. Thus, this confirms that the Cp of Thir can 
be used as a phonon blank to determine the electronic 
contribution to the heat capacity of Npir. 


where Wc is the unit cell volume. The Griineisen param¬ 
eter provides an estimation of the anharmonicity strength 
in a compound. 

The calculated total lattice thermal conductivity of 
Npir is shown in Fig. (TU] From temperatures of 30 K 
to 100 K an exponential decrease of kl{T) is observed, 
which is due to the exponential increase of the phonon- 
phonon scattering via the Umklapp mechanism. For tem¬ 
peratures above 100 K the Umklapp mechanism governs 
the scattering processes and consequently an intrinsically 
low thermal conductivity arises. The lattice thermal con¬ 
ductivity at room temperature assumes a value of 0.64 
Wm“^K“^ and a value of 0.19 Wm“^K“^ at 970 K. 
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Note that these are ultralow values^; for comparison, re¬ 
cent measurements on orthorhombic SnSe crystals with 
a very high thermoelectric figure of merit gave room- 
temperature values between 0.5 and 0.7 Wm“^K“^, and 
values of 0.23 — 0.34 Wm“^K“^ at 970 K, depending on 
the crystal axis^. Very recent ab initio calculations for 
NaBi predicted ultralow values of about 2 Wm“^K“^ at 
300 K^. Npir is thus predicted to have a record low lat¬ 
tice thermal conductivity at high temperatures. Apart 
from a low lattice thermal conductivity, a high electrical 
conductivity is desirable, too, for suitable thermoelec¬ 
tric materials^. As an intermetallic, NpIr is expected 
to have a good electrical conductivity and also a consid¬ 
erably larger contribution to the electronic thermal con¬ 
ductivity than in the chalcogenide systems. However, it 
is not currently possible to simply distinguish between 
the phonon and electron contributions experimentally. 

Low lattice thermal conductivities can be found for 
compounds with a large molecular weight or a complex, 
anisotropic crystal structures^; both conditions are ful¬ 
filled for Npir. In Fig. [in] we in addition show the 
axis-projected thermal conductivities as well as the off- 
diagonal components (in the inset). The latter arise be¬ 
cause of the low symmetry of the monoclinic structure. 
The three crystallographic axis-projected thermal con¬ 
ductivities are of similar size in Npir. 

To assess the importance of the lattice anharmonicities 
for the low thermal conductivity the Griineisen parame¬ 
ters are evaluated. The calculated q-averaged Griineisen 
parameters projected on the crystallographic axes are: 
7 q = 2.46, 7 b = 3.69, and 7 c = 2.46. As the q- 
dependent 7 „q values can be negative, their absolute val¬ 
ues have been computed. The values for Npir are l arg e 
and anisotropic (comparable to those for SnSe, Ref. 14411 , 
which provides evidence for substantial lattice anhar¬ 
monicities that induce heat dissipation and low values 
of the thermal conductivity. In addition, the Griineisen 
parameter 7 b along the 6-axis is much larger than those 
along the a or c-axis. From this we can infer that the 
phonon modes along the 6-axis are more strongly anhar- 
monic and this leads to a weak interatomic bonding and 
hence a good channel to dissipate phonon transport along 
the 6-axis. 


V. CONCLUSIONS 

In conclusion, a new binary equiatomic Np-Ir inter¬ 
metallic has been successfully synthesized. Although it 
is isostructural with Ulr, it is found to be paramag¬ 
netic down to 0.55 K, despite the presence of possibly 
antiferromagnetic interactions. The effective paramag¬ 
netic moment of 3.20 ± 0.02 ^s/Np does not agree well 
with estimates for a free Np^+ or Np^+ ion in either the 
Russell-Saunders or intermediate coupling schemes, im¬ 
plying some degree of 5/ delocalisation; whilst heat ca¬ 
pacity measurements indicate that Npir is a moderate 
heavy fermion system. The form of the electrical resis- 



T(K) 


FIG. 10: (Color online) Calculated lattice thermal conductiv¬ 
ity K of Npir as a function of temperature. Shown is the total 
thermal conductivity kl as well as the thermal conductivities 
along the crystallographic axes. The off-diagonal components 
of the thermal conductivity are given in the inset. 

tivity as a function of temperature and a low temperature 
Sommerfeld coefficient, which is strongly enhanced rela¬ 
tive to that obtained at high temperatures, both indicate 
that Npir should be regarded as a Kondo system. 

Ab initio calculations reveal that the GGA-|-t/ approx¬ 
imation provides a better description of the structural 
and electronic properties of Npir than the plain GGA ap¬ 
proach. After relaxation, the calculations give the same 
geometrical structure as the experimental one; however, 
the calculations without spin-orbit interaction suggest 
that ferromagnetic order is energetically the most favor¬ 
able, followed by the paramagnetic and antiferromagnetic 
ordered states. This result stands in contrast to the lack 
of any experimental observations of ferromagnetism, and 
suggest that DLM calculations with spin-orbit interac¬ 
tion are needed to address this issue thoroughly. The ab¬ 
solute value of the local spin magnetic moments on the 
Np ions is of the order 3.7 /tb. However, due to a size¬ 
able opposite orbital magnetic moment, when including 
the spin-orbit interaction, the moment drops to 2.81 /tb, 
which of a similar magnitude to that extracted from mag¬ 
netic susceptibility measurements. The calculated lattice 
heat capacity of Npir is in good agreement with the mea¬ 
sured heat capacity of Thir, which hence can be regarded 
as a phonon blank for Npir. The lattice thermal conduc¬ 
tivity of Npir is predicted to be exceptionally low at high 
temperatures. 
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Appendix: Optimized structure of Npir 

In Table HIl we give the GGA-I-C/ optimized atomic po¬ 
sitions of the NpIr compound, which was computed to 


crystallize in the monoclinic structure. 


TABLE II: Calculated atomic positions for Npir (space group 
P2i). 


atom 

Wyckoff 

X 

y 

z 

Npl 

2 a 

0.12498 

-0.01371 

0.12513 

Np2 

2 a 

0.62511 

-0.01370 

0.62492 

Np3 

2 a 

0.87507 

0.71392 

0.37498 

Np4 

2 a 

0.37500 

0.71391 

0.87516 

Irl 

2 a 

0.12498 

0.25937 

0.12504 

Ir2 

2 a 

0.62492 

0.25939 

0.62481 

Ir3 

2 a 

0.87489 

0.44092 

0.37484 

Ir4 

2 a 

0.37503 

0.44089 

0.87510 
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